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ABSTRACT 

The structure of e"'"e~ — > bbg events was studied using decays recorded in the SLD 
experiment at SLAG. Three-jet final states were selected and the CCD-based vertex 
detector was used to identify two of the jets as b or b. Distributions of the gluon 
energy and polar angle were measured over the full kinematic range, and compared 
with perturbative QCD predictions. The energy distribution is potentially sensitive to 
an anomalous b chromomagnetic moment k. We measured k to be consistent with zero 
and set limits on its value, —0.11 < n < 0.08 at 95% c.l. (preliminary). 
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The observation of e^e~ annihilation into final states containing three hadronic jets, 
and their interpretation in terms of the process e"^e" qqg provided the first 
direct evidence for the existence of the gluon, the gauge boson of the theory of strong 
interactions, Quantum Chromodynamics (QCD). In subsequent studies the jets were 
usually energy ordered, and the lowest-energy jet was assigned as the gluon; this is 
correct roughly 80% of the time, but preferentially selects low-energy gluons. If the 
gluon jet could be tagged explicitly, event-by-event, the full kinematic range of gluon 
energies could be explored, and more detailed tests of QCD could be performed 0. 
Due to advances in vertex-detection this is now possible using e"^e^ bbg events. The 
large mass and relatively long lifetime, ~ 1.5 ps, of the leading B hadron in 6-quark 
jets 1^ lead to decay signatures which distinguish them from lighter-quark (m, d, s or 
c) and gluon jets. We used our original (1992-5) and upgraded (1996-8) CCD vertex 
detectors 0, ^ to identify in each event the two jets that contain the B hadrons, and 
hence to tag the gluon jet. This allowed us to measure the gluon energy and polar-angle 
distributions over the full kinematic range. 

Additional motivation to study the bbg system has been provided by measurements 
involving inclusive bb decays. Several reported determinations |0 of i?b = 

r{Z^ -^bb) /T{Z^ -^qq) and the Z^-b parity-violating coupling parameter, A^, differed 
from Standard Model (SM) expectations at the few standard deviation level. Since 
one expects new high-mass- scale dynamics to couple to the massive third-generation 
fermions, these measurements aroused considerable interest and speculation. We have 
therefore investigated in detail the strong-interaction dynamics of the 6-quark. We 
have compared the strong coupling of the gluon to 6-quarks with that to light- and 
charm-quarks 0, as well as tested parity (P) and charge©parity (CP) conservation at 
the hbg vertex We have also studied the structure of hbg events via the distributions 
of the gluon energy and polar angle with respect to (w.r.t.) the beamline 0; here we 
present a preliminary update of these measurements using a data sample more than 
3 times larger than in our earlier study. We compare these results with perturbative 
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QCD predictions, including a recent calculation at next-to-leading order (NLO) which 
takes quark mass effects into account [|10[] . 

In QCD the chromomagnetic moment of the b quark is induced at the one-loop level 



and is of order a^/ tt- A more general bbg Lagrangian term with a modified coupling |11 
may be written: 

C''^ = gsbTai^, + '-^{k - tkj,)}bG^^, (1) 

where k and k parameterize the anomalous chromomagnetic and chromoelectric mo- 
ments, respectively, which might arise from physics beyond the SM. The effects of the 
chromoelectric moment are sub-leading w.r.t. those of the chromomagentic moment, 
so for convenience we set k to zero. A non-zero k would modify [|ll|] the gluon energy 
distribution in bbg events relative to the standard QCD case. Hence we have used our 
larger data sample to set improved limits on k. 

We used hadronic decays of bosons produced by e~^e~ annihilations at the SLAC 
Linear Collider (SLC) which were recorded in the SLC Large Detector (SLD) ||12|. The 
criteria for selecting Z° decays, and the charged tracks used for flavor-tagging, are 
described in [^, [l^. We applied the JADE algorithm [|I4| to define jets, using a scaled- 
invariant-mass criterion ycut = 0.02. Events classified as 3-jet states were retained if 
all three jets were well contained within the barrel tracking system, with polar angle 
I cosdjetl < 0.71. From our 1993-98 data samples, comprising roughly 500,000 hadronic 

decays, 126,871 events were selected. In order to improve the energy resolution the 
jet energies were rescaled kinematically according to the angles between the jet axes, 
assuming energy and momentum conservation and massless kinematics [§]. The jets 
were then labelled in order of energy such that Ei > E2 > E^. 

Charged tracks with a large transverse impact parameter w.r.t. the measured in- 
teraction point (IP) were used to tag hbg events 0. The resolution on the impact 



parameter, projected in the plane normal to the beamline, rf, is ad = ll©70/(p_L-\/sin^) 



(1993-5) and 8©33/(pxVsin^) (1996-8) yum, where p± is the track transverse momen- 
tum in GeV/c, and 6 the polar angle, w.r.t. the beamline. The flavour tag was based 
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on the number of tracks per jet, Nlfg, with d/ad > 3 [|. Events were retained in 
which exactly two jets were 6-tagged by requiring each to have > 2, and in which 
the remaining jet had N^^g < 2 and was hence tagged as the gluon; 8196 events were 
selected. The efficiency for selecting true bbg events is 12%. This was estimated using a 



simulated event sample generated with JETSET 7.4 [|15l, with parameter values tuned 



to hadronic e~^e~ annihilation data [jl6[, combined with a simulation of S-decays tuned 
to T(4S) data |]T^ and a simulation of the detector. The efficiency peaks at about 15% 
for 15 GeV gluons. Lower-energy gluon jets are sometimes merged with the parent 
6-jet by the jet-finder. At higher gluon energies the the correspondingly lower-energy 
6-jets are harder to tag, and there is also a higher probability of losing a jet outside 
the detector acceptance. 

For the selected event sample. Fig. 1 shows the N^^g distributions separately for 
jets 1, 2 and 3. In about 16% of cases the gluon-tagged jet is not the lowest-energy 
jet (jet 3). The simulated contributions from true gluons are indicated and the 



estimated gluon purities [|l^] are listed in Table 0. The inclusive gluon purity of the 
tagged-jet sample is 93%. With this sample we formed the distributions of two gluon- 
jet observables, the scaled energy Xg = 2Egi^on/ ^/s, and the polar angle w.r.t. the 
beamline, 6g. The distributions are shown in Fig. 2. The simulation is also shown; it 
reproduces the data. 

The backgrounds were estimated using the simulation and are of three types: non-66 
events, bb but non-bbg events, and true bbg events in which the gluon jet was mis-tagged 
as a 6-jet. These are shown in Fig. 2. The non-66 events (~ 5% of the bbg sample) 
are mainly ccg events, 92% of which had the gluon correctly tagged. There is a small 
contribution (~ 0.1% of the bbg sample) from light-quark events. The dominant back- 
ground is formed by 66 but non-bbg events. These are true 66 events which were not 
classified as 3-jet events at the parton level, but which were mis-reconstructed and 
tagged as 3-jet bbg events in the detector using the same jet algorithm and ycut value. 
This arises from the broadening of the particle flow around the original 6 and 6 di- 
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rections due to hadronisation and the high-transverse-momentum 5-decay products, 
causing the jet-finder to reconstruct a 'fake' third jet, which is almost always assigned 
as the gluon. The population of such fake gluon jets peaks at low energy (Fig. 2(a)), 
as expected. Mis-tagged events comprise less than 1% of the bbg sample. 

The distributions were corrected to obtain the true gluon distributions D*™'^(X) by 
applying a bin-by-bin procedure: L)*™^(X) = C{X) (L)^'''"(X) - S(X)), where X = Xg 
or cos6g, D™'^{X) is the raw distribution, B{X) is the background contribution, and 
C{X) = D^ll^^{X) / D^^f^^i^X) is a correction that accounts for the efficiency for accept- 
ing true hbg events into the tagged sample, as well as for bin-to-bin migrations caused 
by hadronisation, the resolution of the detector, and bias of the jet-tagging technique. 
Here D^^^{X) is the true distribution for MC-generated hbg events, and D\}q'^(X) is 
the resulting distribution after full simulation of the detector and application of the 
same analysis procedure as applied to the data. 

As a cross-check, an alternative correction procedure was employed in which bin- 
to-bin migrations, which can be as large as 20%, were explicitly taken into account: 
L)*™e(^Xi) = M(Xi, Xj) (L)™"'(Xj)-5(Xj))/e(X,), with the unfolding matrix M(Xi, Xj) 
defined by D^7c(^*) = ^i^i^ ^j)DMc"'i-^j)^ where true bbg events generated in bin 
i may, after reconstruction, be accepted into the tagged sample in bin j. e{X) is 
the efficiency for accepting bbg events in bin i into the tagged sample. The result- 
ing distributions of Xg and cos^^ are statistically indistinguishable from the respective 
distributions yielded by the bin-by-bin method. 

The fully-corrected distributions are shown in Fig. 3. Since, in an earlier study ||^, 
we verified that the overall rate of bbg-event production is consistent with QCD expec- 
tations, we normalised the gluon distributions to unit area and we study further the 
distribution shapes. The Xg distribution rises, peaks around Xg ~ 0.15, and decreases 
towards zero as — > 1. The peak is a kinematic artifact of the jet algorithm, which 
ensures that gluon jets are reconstructed with a non-zero energy which depends on the 
yc value. The cos^^ distribution is fiat. 
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We have considered sources of systematic uncertainty that potentially affect our 
results. These may be divided into uncertainties in modelling the detector and un- 
certainties in the underlying physics modelling. To estimate the first case we sys- 
tematically varied the track and event selection requirements, as well as the tracking 
efficiency |0, [T^. In the second case parameters used in our simulation, relating mainly 
to the production and decay of charm and bottom hadrons, were varied within their 
measurement errors [0. For each variation the data were recorrected to derive new Xg 
and cos^g distributions, and the deviation w.r.t. the standard case was assigned as a 
systematic uncertainty. None of the variations affects our conclusions. All uncertain- 
ties were conservatively assumed to be uncorrelated and were added in quadrature in 
each bin of Xg and cos6g. 

We compared the data with perturbative QCD predictions for the same jet algo- 
rithm and i/c value. We used leading-order (LO) and NLO results based on recent 
calculations |]lO| in which quark mass effects were explicitly taken into account; a b- 
mass value of mi,{mz) = 3GeV/c^ was used [|^. We also derived these distributions 
using the 'parton shower' (PS) implemented in JETSET. This is equivalent to a calcu- 
lation in which all leading, and a subset of next-to-leading, Ini/c terms are resummed 
to all orders in as- In physical terms this allows events to be generated with multiple 
orders of parton radiation, in contrast to the maximum number of 3 (4) partons al- 
lowed in the LO (NLO) calculations, respectively. Configurations with > 3 partons are 
relevant to the observables considered here since they may be resolved as 3-jet events 
by the jet-finding algorithm. 

These predictions are shown in Fig. 3. The calculations reproduce the measured 
cos6g distribution, which is clearly insensitive to the details of higher-order soft parton 
emission. For Xg, although the LO calculation reproduces the main features of the 
shape of the distribution, it yields too few events in the region 0.2 < Xg < 0.5, and too 
many events for < 0.1 and Xg > 0.5. The NLO calculation is noticeably better, but 
also shows a deficit for 0.2 < < 0.4. The PS calculation describes the data across the 
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full Xg range. These results suggest that multiple orders of parton radiation need to be 
included, in agreement with our earlier measurements of jet energy distributions using 
flavor- inclusive decays pO| . We also investigated LO and NLO predictions based on 
matrix elements implemented in JETSET which assume massless quarks. The resulting 
distributions are practically indistinguishable from the massive ones, even though the 
large 6-mass has been seen [|19| to affect the bbg event rate at the level of 5%. The 
effect of varying as within the world-average range is similarly small. 

We conclude that perturbative QCD in the PS approximation accurately repro- 
duces the gluon distributions in bbg events. However, it is interesting to consider the 
extent to which anomalous chromomagnetic contributions are allowed. The Lagrangian 
represented by Eq. □ yields a model that is non-renormalisable. Nevertheless tree-level 



predictions can be derived |Tl|] and used for a 'straw man' comparison with QCD. 
For illustration, the effect of a large anomalous moment, k = 0.75, on the shape of 
the Xg distribution is shown in Fig. 3(a); there is a clear depletion of events in the 
region Xg < 0.5 and a corresponding enhancement for Xg > 0.5. By contrast the shape 
of the cosOg distribution is relatively unchanged (not shown), even by such a large k 
value. In each bin of the Xg distribution, we parametrised the leading-order effect of 
an anomalous chromomagnetic moment and added it to the PS calculation to arrive 
at an effective QCD prediction including the anomalous moment at leading-order. A 
minimisation fit was performed to the data with k as a free parameter, yielding 
K, = —0.011 ± 0.048(stat.)]J^Q QQ3(syst.), which is consistent with zero within the errors, 
with a of 17.8 for 9 degrees of freedom. The distribution corresponding to this fit 
is indistinguishable from the PS prediction (Fig. 3(a)) and is not shown. Our result 
corresponds to 95% confidence- level (c.l.) upper limits of —0.11 < k < 0.08 (prelimi- 
nary) . 

In conclusion, we used the precise SLD tracking system to tag the gluon in 3- 
jet e^e^ ^ ^ bbg events. We studied the structure of these events in terms of 
the scaled gluon energy and polar angle, measured across the full kinematic range. We 
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compared our data with perturbative QCD predictions, and found that the effect of the 
6-mass on the shapes of the distributions is small, that beyond-LO QCD contributions 
are needed to describe the energy distribution, and that the parton shower prediction 
agrees best with the data. We also investigated an anomalous 6-quark chromomagnetic 
moment, k, which would affect the shape of the energy distribution. We set preliminary 
95% c.l. limits of -0.11 <k< 0.08. 
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Jet label 


# Tagged gluon jets 


Purity 


3 


1140 


94.4 % 


2 


155 


90.1 % 


1 


34 


73.1 % 



Table 1: Estimated purities of the tagged gluon-jet samples. 



QCD Calculation 


X^: a^g (10 bins) 


LO rribimz) = 3 GeV/c^ 


170 


NLO mb(mz) = 3 GeV/c^ 


51 


PS Mb = 5 GeV/c^ 


18 



Table 2: fo^^ the comparison of the QCD predictions with the corrected data. 
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Figure 1: The N^^g distributions for jets in bbg-t&gged events, labelled according to jet 
energy (dots); errors are statistical. Histograms: simulated distributions showing jet 
flavour contributions. 
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Figure 2: Raw measured distributions of (a) Xg and (b) cos^^ (dots); errors are statis- 
tical. Histograms: simulated distributions including background contributions. 
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Figure 3: Corrected distributions of (a) Xg and (b) cos^^ (dots); errors are statistical. 
Perturbative QCD predictions (see text) are shown as lines joining entries plotted at 
the respective bin centers. 
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